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FURTHER EXPERIMENTAL RESULTS ON THE INFLUENCE OF THE 


TURBULENT BOUNDARY LAYER ON PANEL FLUTTER 


Peter A. Gaspers, Jr., Lado Muhlstein, Jr., 
and Daniel N. Petroff 

Ames Research Center 


SUMMARY 


Flutter boundaries were obtained for a flat rectangular panel for ratios 
of boundary- layer thickness to panel length of 0.036 to 0.128 at Mach numbers 
of 1.2 to 1.4. The panel was isotropic and unstressed with a length-to-width 
ratio of 0.5 and the edges were rotational ly restrained to approximate a 
clamped condition. The results show that the turbulent boundary layer has a 
large stabilizing influence. The effect is most pronounced at a Mach number 
of 1.2 and decreases rapidly as the Mach number increases. The experimental 
results, extrapolated to zero boundary- layer thickness, are compared with the 
available theoretical predictions. A theoretical analysis of the vibration of 
plates with variable rotational edge restraint is given and, within its 
framework, the edge restraint parameters for the experimental panels are 
determined. 


INTRODUCTION 


Experimental data in a previous report (ref. 1) showed that the 
turbulent boundary layer has a large stabilizing influence on panel flutter at 
low supersonic Mach numbers. In that study the Mach numbers ranged from 1.05 
to 1.4, the stiffness parameter was about 0.072, and the mass ratio parameter 
ranged from 0.06 to 0.11. In the present study additional experimental 
results have been obtained for Mach numbers from 1.2 to 1.4 but for a stiff- 
ness parameter of about 0.025 and mass ratios of 0.009 to 0.017. The panel 
made of Invar was flat, rectangular (length-to-width ratio 0.5), and 
unstressed with rotationally restrained, approximately clamped, edges. Exten- 
sive precautions, discussed in detail in reference 1, were taken to eliminate 
as many as possible of the usual sources of experimental inaccuracy. The 
results from reference 1 and from the present investigation, including extrap- 
olations to zero boundary- layer thickness, are compared with theoretical pre- 
dictions based on three-dimensional unsteady potential flow. The present 
results support the findings of reference 1 and show that the influence of the 
boundary layer can account for a substantial portion of the disagreement 
between theoretical and experimental panel flutter results in the low 
supersonic Mach number range (refs. 1-4). 



APPARATUS 


Panel Configuration 

Dimensions and mounting details of the panel and frame used in this 
experiment are shown in figure 1. The 9.00-inch by 18.00-inch by 0.0192-inch- 
thick panel was made from Invar, an iron-nickel alloy (64-percent iron, 36- 
percent nickel) with a very low coefficient of thermal expansion. The panel 
was attached to an Invar frame, 11.00 inch by 20.00 inch by 0.340 inch thick, 
with an adhesive (Eastman 910) . 

Closely spaced roundhead screws were threaded into the panel sheet from 
the cavity side (air stream side smooth) to provide additional clamping action 
and to prevent loss of the panel in case of bond failure during flutter. Com- 
parison of panel natural frequencies before and after the flutter test 
disclosed no deterioration of the adhesive bond. 

The panel frame was attached to the test fixture by one high-stiffness 
and two low-stiffness flexures as shown in figure 1. This three-point mount- 
ing was designed to minimize the introduction of external stresses into the 
panel frame. 


Wind Tunnel 

The experiments were conducted in the Ames 2- by 2-Foot Transonic Wind 
Tunnel, which is of the continuous flow type with porous test section walls 
surrounded by a plenum chamber (see ref. 5) . Mach number is continuously 
variable from 0.60 to 1.40 and dynamic pressure is continuously variable from 
200 to 1500 psf. 


Variable Boundary- Layer Test Fixture 

The fixture for varying the boundary- layer thickness consisted of a 
splitter plate with a sharp leading edge installed in one side wall of the 
wind tunnel. This splitter plate could be positioned by remote control to be 
flush with the wall or to project up to 1.0 inch into the airstream. Dimen- 
sions of the test fixture are shown in figure 2. When the splitter plate was 
flush with the tunnel wall, the boundary- layer thickness in the test region 
was maximum. When the splitter plate was moved into the airstream, a portion 
of the tunnel-wall boundary layer was directed into the plenum chamber sur- 
rounding the test section and a thinner boundary layer was established on the 
splitter plate. The boundary- layer thickness in the panel test region is var- 
iable by this means between the approximate limits of 0.25 and 1.00 inch. 
Boundary- layer thickness was measured by four retractable boundary-layer 
total -pressure probes located as shown in figure 2. Figures 3 and 4 are front 
and rear views and show tunnel installation of the variable boundary- layer 
test fixture. 
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Other Equipment 


The panel was backed by a sealed rectangular cavity (identical to that 
used in ref. 1) with the same dimensions as the panel and with an effective 
depth of 7.5 inches (where effective depth is defined as the actual cavity 
volume minus volume occupied by instrumentation all divided by panel area) . 

The panel frame was sealed to the cavity with adhesive plastic tape to minimize 
the possibility of introducing stresses into the panel frame. 

The differential pressure between the cavity and free stream was 
controlled manually to within ±0.25 psf by a system similar to that described 
in reference 2. 


Panel Natural Frequencies, Material Properties, 
and Boundary Conditions 

In order to calculate the natural frequencies of a plate accurately, the 
plate dimensions, elastic constants, density, and boundary conditions must be 
known accurately. The plate dimensions and density can be measured to almost 
an arbitrary degree of accuracy and so present no problem. Determining the 
elastic constants and boundary conditions, however, presents difficulties 
since the two effects are not easily separated. One method (ref. 6) of deter- 
mining the elastic constants of thin isotropic plates depends on the fact that 
the normalized frequency spectrum of a plate with free edges is a function of 
Poisson’s ratio. Therefore, for the free plate the experimentally obtained 
ratio of any two natural frequencies together with the value of any natural 
frequency determines Poisson’s ratio and Young’s modulus. The free circular 
plate is most convenient for this purpose because its eigenvalues can be 
calculated exactly. 

This method was used here to redetermine the elastic constants of the 
magnesium alloy sheet material used in reference 1 and the results are shown in 
table 1. The resulting Young’s modulus of 6.2xl0 6 psi lies between our pre- 
vious measured value of 5.6x10 s psi and the handbook value of 6.5x10 s psi. 

The value of Poisson’s ratio coincides with the handbook value. When the same 
method was applied to the Invar sheet material, an unacceptably large amount 
of data scatter was obtained and so, for that case, the handbook values are 
used. The failure of this method for the Invar sheet material may be due to 
residual stresses or small variations in thickness or flatness. The eigen- 
values of a free plate depend rather weakly on Poisson’s ratio and any varia- 
tions between the experimental specimen and the theoretical model can be 
expected to cause difficulty. 

The boundary conditions of the magnesium panel and of the Invar panel 
were determined within the framework of an analytical model that prescribed 
zero transverse edge displacement but arbitrary rotational edge restraint. A 
theoretical analysis of this model together with some typical numerical results 
is presented in the appendix. The same theoretical model is analyzed in ref- 
erence 7 by a different method. It was assumed that =q^ =q^ and 

c l r)0 = q - (equal rotational restraints on opposite panel edges); then 

q^ and q^ were determined to give a least squares fit to the experimental 
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natural frequencies. The resulting values of q ^ and q n (see table 1) show 
that the magnesium panel was considerably closer to the clamped condition than 
the Invar panel. The lower values of the restraint parameter for the Invar 
panel are probably associated with the adhesive bonded joint between panel and 
frame (the magnesium panel of ref. 1 was an integral part of the frame) , but 
the detailed mechanism of declamping is not clear. The ratio of q n to q^ 
for both panels is much larger than the theoretical value of 2 which results 
if the dimensional rotational spring constants 0x and 9y on all four sides 
are equal. This may be due to the variation of 6 X and 0y along the edge 
caused by the torsional and bending stiffness of the frame. It should be 
noted, however, that as both q^ and q n become large, the panel approaches 
the clamped condition and the ratio of q^ to q^ becomes less significant. 

The values of q^ and qn given in table 1 result in a root-mean-square 
error between theoretical and experimental natural frequencies (the first 5) 
of 1.3 percent for the magnesium panel and 2.4 percent for the Invar panel. 

If clamped edges are assumed the RMS errors become 3.6 percent for the 
magnesium panel and 11.9 percent for the Invar panel. 


Thermal Effects 

As discussed in some detail in reference 1, a relatively small 
temperature difference between the panel and frame can have a large effect on 
the natural frequencies and consequently on the flutter boundary of thin 
panels. To prevent excessive temperature sensitivity of the low stiffness 
panel desired for this test, it was constructed of Invar which has a coeffi- 
cient of thermal expansion approximately l/20th that of aluminum or magnesium. 

The theoretical temperature differential (assuming a clamped panel) 
between the panel and frame that would buckle this panel is 17.5° F. An 
experiment described in reference 1 was unable to produce buckling and the 
change in natural frequency observed was only a fraction of a hertz per °F. 

The maximum temperature differential between panel and frame was limited to 
0.5° F during acquisition of flutter data. Therefore, the influence of 
temperature differential can be assumed negligible. 


INSTRUMENTATION 


The instruments used and measurements made during this test are identical 
to those reported in reference 1. The panel displacement was measured with a 
noncontacting capacitance- type transducer located as shown in figure 2. The 
surface static pressures and the total pressures from the boundary- layer 
probes were measured with conventional strain-gage differential-pressure 
transducers connected to the orifices by short lengths of tubing and refer- 
enced to free-stream static pressure. Panel and frame temperatures were mea- 
sured with iron-constantan thermocouples in the center of the panel and at two 
points on the frame. The thermocouple in the center of the panel was con- 
structed of wire 0.005 inch in diameter and was the only object contacting the 
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panel. Steady quantities such as temperatures and pressures were recorded 
digitally. Unsteady or slowly varing quantities such as fluctuating static 
pressure, panel response, and pressure differential between the cavity and 
free stream were recorded by a frequency modulated magnetic tape recorder with 
a flat frequency response (±0.5 dB) from dc to 1250 Hz. 


TEST PROCEDURE 

Flow Conditions in Panel Test Region 

The flow conditions in the panel test region are an important 
consideration in panel flutter testing, as discussed in reference 1. Both the 
surface static-pressure distribution and the fluctuating pressure were mea- 
sured. The maximum variation of surface static-pressure coefficient (Cp) is 
±0.02, with most variations being less than ±0.01. The power spectra of the 
fluctuating pressures indicate that there are no dominant frequencies in the 
range of interest for this investigation. Plots of typical static-pressure 
distributions, fluctuating pressure coefficients, power spectra of the fluctu- 
ating pressures, and boundary- layer velocity profiles for various splitter 
plate displacements are presented in reference 1. 


Determination of the Flutter Boundary 

The many problems involved in accurately defining an experimental flutter 
boundary are discussed in detail in reference 1. In the present study, the 
flutter boundary is determined in the same manner as in reference 1, namely, 
which is as follows: At dynamic pressures well below that for flutter, the 

panel response to pressure disturbances from the turbulent boundary layer 
increases approximately linearly with dynamic pressure, as indicated in fig- 
ure 5. In this figure, panel amplitude in the mode which becomes the flutter 
mode is shown as a function of dynamic pressure. This ’’subflutter response” 
is used as a reference level to define the flutter point systematically. The 
straight lines labeled 2x, 3x, etc., in figure 5 are multiples of the subflut- 
ter response. Their intersections with the amplitude response curve define 
values of dynamic pressure necessary to produce panel amplitudes which are 
these multiples of the subflutter response. Figure 6 shows dynamic pressure 
as a function of the multiples of the subflutter response for the same condi- 
tions as shown in figure 5. As the multiple increases, the difference in 
dynamic pressure between successive multiples decreases; that is, the dynamic 
pressure appears to be converging to a limit. Furthermore it appears that the 
6x multiple results in a dynamic pressure very close to this limiting dynamic 
pressure and is therefore arbitrarily defined herein as the flutter point. 

The following test procedure was adhered to throughout the program. Mach 
number and splitter plate position were held constant while dynamic pressure 
was increased in small increments beginning at a value well below the flutter 
boundary. After each increase in dynamic pressure, the panel-frame tempera- 
ture differential was allowed to stabilize to ±0.5° F to minimize thermal 
stress effects. The indicated differential pressure between the cavity and 
free stream was then varied slowly from -50 to +50 psf while the panel 
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response was recorded. It is assumed that the maximum response so measured 
corresponds to zero average differential pressure across the panel. Next the 
indicated differential pressure was set to the value at which maximum response 
occurred. This response was recorded so that amplitude spectra could be 
obtained with all test conditions fixed. The above procedure was repeated for 
each increase in dynamic pressure until the maximum response was a large- 
amplitude, essentially sinusoidal motion indicating that the flutter boundary 
had been penetrated. Boundary- layer velocity profiles were obtained at 
dynamic pressures near the flutter boundary. The forward probe was retracted 
while data were being recorded from the rear probes, and all probes were 
retracted when boundary- layer data were not being recorded. 

The maximum response amplitude of the panel mode which becomes the flutter 
mode was determined from plots of response amplitude versus differential pres- 
sure for each dynamic pressure with the special data reduction system described 
in reference 1. 


RESULTS AND DISCUSSION 


Influence of the Turbulent Boundary Layer 
on Flutter Dynamic Pressure 

The dynamic pressures for various multiples of subflutter response for a 
series of boundary- layer thicknesses are shown in figure 7 for Mach numbers of 
1.20, 1.30, and 1.40. The geometric boundary- layer thickness used is the 
thickness in the center of the panel, estimated by assuming linear growth 
between the upstream and downstream probes. The curve for 6x subflutter 
response is considered the flutter boundary. The flutter dynamic pressure in 
figure 7 is approximately a linear function of the geometric boundary- layer 
thickness, thus allowing reasonable extrapolation to zero boundary- layer 
thickness . 

In figure 8 the flutter dynamic pressure normalized by the extrapolated 
flutter dynamic pressure for zero boundary- layer thickness at that Mach number 
(q/q<5 = o) pl° t1:e d as a function of the ratio of boundary- layer thickness to 
panel length (6/a). The slopes of the curves in figure 8 d(q/q£ =0 )/d(6/a) 

are plotted as a function of Mach number in figure 9, which shows the strong 
dependence of boundary- layer effects on Mach number. The maximum effect 
occurs below M = 1.20, the minimum Mach number in this experiment, and is 
similar to the result reported in reference 1. 


Influence of the Turbulent Boundary Layer on Flutter Frequency 

The flutter frequency normalized by the flutter frequency for zero 
boundary- layer thickness (f/f^o)* obtained by linear extrapolation in the 
manner of reference 1, is plotted versus the normalized boundary- layer thick- 
ness in figure 10. Note that an increase in boundary- layer thickness causes a 
decrease in flutter frequency. 
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The slopes of the curves in figure 10, d(f/f,5 = o)/d(6/a) > are presented 
in figure 11 as a function of Mach number and indicate that the minimum effect 
of boundary- layer thickness on flutter frequency for this panel occurs near 
M = 1.30. 

The results shown in figure 10 are replotted in figure 12 where the 
reduced frequency of flutter is presented as a function of Mach number for 
various normalized boundary- layer thicknesses. 


Generalized Panel Flutter Boundaries 

A natural set of dimensionless parameters for the flutter of flat 
isotropic panels exposed to linearized three-dimensional unsteady potential 
flow (ref. 8) is the mass ratio parameter, y = p^a/p^h, the stiffness param- 
eter, k = /r 2 E/C OQ 2 p s (1 - v 2 ), the panel length-to-width ratio, a/b, the Mach 
number M, and the reduced frequency k . Flutter boundaries, along which the 
reduced frequency varies, may be plotted in the k-y plane for each combina- 
tion of Mach number and length-to-width ratio. Various authors have replaced 
k and y with other sets of dimensionless parameters (see refs. 2, 9, 10) 
but all of these alternate sets depend on Mach number, length-to-width ratio, 
or panel edge conditions and so, as pointed out in reference 8, plots for dif- 
ferent Mach numbers, etc., cannot be directly compared. For this reason all 
flutter boundaries presented herein are plotted in the k-y plane. If 
boundary- layer effects are to be considered, the ratio of boundary- layer 
thickness to panel length, 6/a, and the ooundary- layer profile must also be 
considered. In figure 13, the flutter mass ratio parameter y is plotted as 
a function of boundary- layer thickness to panel length ratio 6/a, for the 
three test Mach numbers. Since these data fall approximately on a straight 
line, they can reasonably be extrapolated to zero boundary- layer thickness 
(6/a = 0) to give a value of y for comparison with results predicted by 
potential flow theory. 


Comparison With Theoretical Predictions 

For the panel length-to-width ratio (a/b) and Mach number range of this 
test, valid panel flutter predictions can be obtained only by three-dimensional 
unsteady aerodynamic theory. Theoretical results based on such theory for 
clamped rectangular panels for the length-to-width ratio and Mach number range 
of interest are presented in references 2 and 10. These results are used for 
comparison with the experimental data presented in this report and in 
reference 1. 

In reference 2, a four-mode analysis was performed for a clamped 
rectangular panel with a/b = 0.46. Results are presented for M = 1.10, 1.20, 
1.30, and /2 and structural damping coefficients of 0 and 0.05 in the dynamic 
pressure parameter-density ratio (X,y) plane where the relation between X, k, 
and y is given by 


X = 12M 2 y/Bk 2 
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It should be noted that the adequacy of the four-mode analysis was based on 
convergence studies using quasi-steady aerodynamics. No convergence studies 
were performed using three-dimensional unsteady aerodynamics. Although these 
results are for a/b = 0.46, they are expected to differ very little from 
those for a/b = 0.5 and thus are used for comparison without correction. In 
another approximation in reference 2 the off-diagonal elements of the stiff- 
ness matrix are neglected, but in reference 11 it is shown that for length-to- 
width ratios in the vicinity of 0.5, the effect on the flutter boundary is 
negligible. For use in this comparison, the A,y plots of reference 2 were 
converted to the k,y plane. A theoretical boundary for M = 1.40 was 
obtained by plotting y versus M for various values of k and obtaining the 
desired y,k plot by interpolation. Since M = 1.40 is not significantly 
different from M = /2, the errors resulting from interpolation are small. 

Reference 10 presents flutter boundaries for clamped panels of a/b = 0.5 
for M = 1.30 only. As in reference 2, the off-diagonal elements of the 
stiffness matrix were neglected. The parameters employed are y^ = 1/y, where 
u is the usual mass ratio parameter (however, 1/yi is plotted), and a stiff- 
ness parameter 


kl 



k 


where a>i and Xi are the lowest natural frequency and eigenvalue, respec- 
tively, of the plate. Results of a six-mode analysis are presented for coef- 
ficients of structural damping of 0 and 0.01 in the ki - 1/yi plane. These 
results have been converted to the k-y plane for use herein. The value of 
\ 1 = 604.05835 used in the conversions was computed by Galerkin's method using 
14 clamped beam modes in both the chordwise and spanwise directions. 

In the following comparisons of the experimental results from reference 1 
and the present study with the theoretical predictions of references 2 and 10, 
the Young's modulus used in the stiffness parameter k would result in a 
first-mode natural frequency for a clamped panel that matches the experimen- 
tally determined first natural frequency. Since the theoretical results are 
for perfectly clamped panels, this choice of E compensates, at least partly, 
for the less than clamped edge conditions of the experimental panels. 

In addition, the stiffness parameter k*, based on the handbook value of 
E, is included in table 2 where the various panel flutter parameters, 
boundary- layer thicknesses, frequencies, and other pertinent quantities are 
tabulated. 

In figure 14(a) the experimental data for M = 1.3 in the k-y plane 
from reference 1 and the present study are compared with theoretical stability 
boundaries from references 2 and 10. The stability boundaries from reference 
10 are more accurate than those of reference 2 since a six-mode analysis was 
used in the former case whereas four modes were used in the latter. Theory 
predicts the existence of several stability boundaries - that giving the 
largest value of k, for a given y, being the critical or flutter boundary. 
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For y corresponding to the experimental results, the theoretical flutter 
boundary is that labeled mode 2 which is predominately composed of the second 
natural, in vacuo, mode. However, the experimental points corresponding to 
zero boundary- layer thickness for both sets of experimental data fall between 
the reference 10 mode 1 boundaries for g = 0 and 0.01. The mode 1 boundary 
is predominately composed of the first natural, in vacuo, mode but contains a 
significant amount of the second natural mode and so it is sometimes called a 
coupled mode instability. The mode 2 boundary is associated with a very weak 
instability that completely disappears for a value of g of about 0.02 (ref. 
2) and was not observed in the present study or in reference 1. The values of 
g for the experimental cases were not measured, but, because of the method of 
panel fabrication, are expected to be in the neighborhood of 0.01 or less. In 
addition, we note from qualitative observations that the damping of the Invar 
panel of the present study is considerably less than that of the magnesium 
panel of reference 1. This is to be expected from thermoelastic considera- 
tions (ref. 12) since the thermal expansion coefficient of Invar is very low. 

The agreement between experimental and theoretical reduced frequencies 
for the mode 1 boundary is considered to be only fair. Since the theoretical 
reduced frequencies associated with the mode 2 boundary are much higher, it is 
clear that the observed flutter is first mode. 

In figure 14(b) the experimental data from reference 1 at M = 1.1 is 
compared to the first-mode theoretical k-y boundaries from reference 2. As 
at M = 1.3, the agreement between theory and experiment improves as the 
boundary- layer thickness approaches zero. However, the difference between 
theoretical and experimental k-y points for zero boundary- layer thickness is 
considerably larger than at M = 1.3. On the other hand, the experimental and 
theoretical reduced frequencies are in good agreement. 

In figure 14(c) the experimental data from reference 1 and the present 
study, for M = 1.2, are compared with the first-mode theoretical boundary 
from reference 2. The effect of the boundary layer is very similar to that 
at M = 1 . 1 and 1.3. 

In figure 14(d) the experimental data at M = 1.4 are compared with the 
theoretical boundaries from reference 2. However, in contrast to the lower 
Mach numbers, as the boundary- layer thickness decreases, the difference 
between theory and experiment increases for both sets of experimental data. 
This behavior is difficult to explain since both the boundary layer and 
structural damping have less effect at M = 1.4 than at lower Mach numbers. 

On the other hand, except at M = 1.3, the agreement between theory and 
experiment for zero boundary- layer thickness is better at M = 1.4 than at the 
lower Mach numbers. 

It is interesting to note that the disagreement between theory and 
experiment is similar for both panels at all Mach numbers, indicating that the 
disagreement results from some systematic inaccuracy or omission in either 
theory or experiment rather than random errors in the experimental results. 

It is possible that a combination of errors, usually considered small individ- 
ually, such as unconverged solutions, neglect of cavity, omission of off- 
diagonal terms in stiffness matrix and imprecise matching of boundary 
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conditions could, taken together, cause an appreciable error in some of the 
theoretical predictions. It should also be remembered that we have used an 
arbitrary, though systematic, definition of the flutter boundary and this 
could cause an error, depending on the strength of the instability. 

Finally, it should be noted that Dowell (ref. 13) previously compared the 
experimental results of reference 1 with the theory from reference 2 and con- 
cluded that experimental data from reference 1 for zero boundary- layer thick- 
ness agree well with the theory. His conclusion is somewhat stronger than 
those made herein and it bears repeating that detailed comparisons at these 
Mach numbers depend strongly on the structural damping. 


CONCLUSIONS 


For the configuration tested the following conclusions can be noted. 

1. The turbulent boundary layer has a large stabilizing influence on 
the flutter of flat isotropic panels at low supersonic Mach numbers. 

2. The effect of the turbulent boundary layer on flutter dynamic 
pressure decreases rapidly with increasing Mach number above M = 1.20 for 
the configuration tested. 

3. The agreement between the experimental results and theoretical 
predictions based on unsteady, three-dimensional, potential flow is consider- 
ably improved when the experimental results are extrapolated to zero boundary- 
layer thickness except at M = 1.4 where the effect of the boundary layer is 
much smaller than at the lower Mach numbers. 

4. The experimental results indicate that flutter occurs as a 
predominantly first-mode instability rather than the predominantly second-mode 
instability predicted by theory. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Jan. 15, 1970 
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APPENDIX A 


NATURAL FREQUENCIES OF RECTANGULAR PLATES WITH 
ROTATIONALLY RESTRAINED EDGES 


The differential equation and boundary conditions for a rectangular 
isotropic plate with arbitrary rotational edge restraint and zero edge 
translation are as follows: 


D(w + 2w + w ) + mw„. = 0 
xxxx xxyy yyyy tt 


Dw 

XX 

- 

0 _w = 0 , 

XO X 

W = 

0 

on 

x = 0 


Dw 

XX 

+ 

0 w = 0 , 
xa x ’ 

w = 

0 

on 

X 

II 


Dw 

yy 

- 

0 w = 0 , 

yo y 

w = 

0 

on 

y = 0 


Dw 

yy 

+ 

0 . w = 0 , 
yb y 

w = 

0 

on 

& 

II 

X 

J 


(Al) 


(A2) 


where 0 xQ , etc., are the rotational spring constants per unit length along 

the edges of the plate and the subscripts on w indicate partial 
differentiation with respect to the subscript. 

Making the transformation: 

x = a 5 y = bn w=aW t=tx 

o 

we obtain the following dimensionless form of the equations: 


W 


5555 


+ 2 

— w__ + — 

b 2 b 4 

w + w 

nnnn 

TT 

LW + W. 

W 

55 

1 

& 

vPY 

O 

vm 

II 

o 

W = 0 

on 

o 

II 

UP 

W 

55 

* '’a”? " 0 • 

W = 0 

on 

ll 

h- 1 

W 

nn 

- %o W n - 0 ' 

W = 0 

on 

o 

n 

w 

nn 

+ %i W n " 0 ' 

W = 0 

on 

n = 1 


TT 


= 0 


-\ 


J 


(A3) 


(A4) 


where 
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a0 


xo 


aG 


xa 


be 


yo 


be 


yb 


4 ?0 


l £l 


l no 


l ni 


are the dimensionless rotational edge restraint parameters and t Q 2 = ma 4 /D. 
Substituting = U(5,r))g(T) into equation (A3) we obtain: 


LU _ fr 
U = " g 


Since the left side is a function only of £ and n while the right side is a 
function of x, they must both be equal to the same constant A and the 
eigenvalue problem is: 


LU - AU = 0 


" q £0 U 5 

= 0 , 

U = 0 

on 

C = 0 

'55 + q U U 5 

= 0 , 

U = 0 

on 

5 = 1 

- q U 

nn ^no n 

= 0 , 

a 

II 

o 

on 

23 

II 

o 

+ q U 
nn ni n 

= 0 , 

U = 0 

on 

n = 1 


g TT + Ag = 0 


(A5) 


(A6) 


(A7) 


From equation (A7) the circular frequencies are given by = (l/t 0 ) An 

approximate solution of the eigenvalue problem can be obtained by Galerkin's 
method as follows: 

Assume a solution in the form: 


J 

U(S,n) = C-jU-jU, n) 

j=i 


where each Uj satisfies the boundary conditions. 
Substituting in equation (A5) we obtain: 
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J 

T: CLUj - AUj)Cj = 0 

j=l 


(A8) 


II 


r 


Now multiply by and integrate over the panel to obtain: 


E / / 

j = l 0 0 


XUjUj )d£ dn Cj = o , 


i = 1, 2, 


. U 


(A9) 


or, in matrix form: 


where 


Let 


where 


(A - XB)C = 0 


(A10) 




a UiLUj d? dn 


0 0 


I 


>ij - S' f 


dc dn 


o o 


(All) 


Uj (5,n) = v m (5)* n (n) 
Ui(5,n) = Mf r (5)+ s Cn) 


(A12) 


%(£) = sin(a^5 + 3^)sin mir? 
<f> n (n) = sin (a^n + 3^) sin nirn 


(A13) 


¥ m and <j> n satisfy the displacement boundary conditions and a ? , 3^, a n , 3^ 
are chosen to satisfy the remaining boundary conditions. 


Substituting equations (A12) and (A13) into (A6) we find that ag, 3^ 
must satisfy the following equations: 


2a^ cos 3^ - q^ 0 sin 3^ = 0 

(A14) 


(2a^ cos + q^ sin a^)cos 3^ + (q^ cos - 2a^ sin a^)sin 3^ = 0 
with corresponding equations for a^, 3^. 


(A15) 
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I 


Eliminating B results in: 

( 4a £ 2 ~ c l ?0 q 5i )tan a ? " 2a 5 (q 50 + q 5i^ = 0 <> 16 ) 

Equation (A16) has one nontrivial root in the interval 0 to it (note that 
ct£ = 0 gives a function that satisfies simple support boundary conditions 
(q = q^ 1 = 0), and = it satisfies clamped conditions (q^ Q = q^ = “)), 

which may be calculated numerically for arbitrary values of q ?Qj q^. 

From equation (A14) , tan B ? = 2a ? /q 50 which determines Bg once is 

known . 

Substituting equations (A12) into (All) we have 


a. . 


il 




t> d> dn + 2 — 

5 n b 2 


d5 f\ <i> lx dn 

J r m J Y s n 


+ 




T r m 



|> (J) 1V 
K s Y n 


dn 



f V m d ? J" dn 

0 0 


The integrals are straightforward, resulting in 


f 

•'o 

f 


Vm d5 = 2D rm + T 6 


rm T 4 u rm 


d £ = (2a r 2 - m 2 u 2 - r 2 u 2 )D rm - j (a^ 2 + m 2 ir 2 )6 m 


/ « = 2 [‘ : 


3a ^ + r^m^TT^ + a 2 fr 2 + m 2 )7r z |D 


)tt 2 ]i 


r rm 


+ j £(a ? + mu) 4 + (a ? - mu) ^6^ 
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where 


0 if r ^ m 


J rm 


1 if r = m 


is the Kronecker delta and 


rmir 2 a^ J^(-l) r+m sin 2 (a^ + B^) - sin 2B^J 


'rm 


16 + (r + m) y - (r + m) j + (r - m) j otg - (r - m) j 


with corresponding integrals for 4> m . At 

<*£ = ±(r ± m) j 

Drm contains removable singularities. If is too close to a singularity, 

numerical accuracy will be impaired in evaluating D rm ; however, in that case 
the singularity may be removed analytically by expanding the numerator of 
Drm in a Taylor series about ± (r ± m)(Tr/2). The result will be: 


Drm - 


16 a 


rimr' 


+ (r + m) 


“{ 2 JT sin ( 26 S * ¥)[“« 1 tr * "° f] 


IT 

2j[ a 5 “ 


(r + m) 


+ (r - m) 


J [ a 5 - Cr - m) J 


The matrix eigenvalue problem (A - XB)C = 0 (eq. (A10)) can be placed in 
standard form in a numerically stable way utilizing a method given in refer- 
ence 14. Factor B into the product TT^ of a lower triangular matrix T 
and its transpose T^. Let C = (tT) -1 Y and we have (A - XTTT) (tT) - 1 Y. Mul- 
tiply by T _1 to give [T - 1 A(tT) -1 - XI]Y = 0 and the problem is in standard 
form. 


NUMERICAL CALCULATIONS 


Eigenvalues and eigenvectors have been calculated by programs available 
from SHARE, NUEIG 5, and GIVHO (refs. 15, 16). In the computations, 10 modes 
were used in the long direction of the panel and 5 modes in the short 
direction. 

In figure 15 /x^ is plotted versus = qr 0 = qgi for a/b = 0.5 and 

q n = 2q^ to indicate the effect of varying the edge condition from simply 


15 



supported to clamped. Note that q n = 2qg corresponds to the case of equal 
dimensional spring constants 9 X = 0y for a/b = 0.5. 

The eigenvalues for perfectly clamped edges as shown in figure 15 were 
calculated by Galerkin's method using 14 clamped beam modes in both the 
spanwise and chordwise directions. 
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TABLE 1.- PANEL NATURAL FREQUENCIES, LOCATIONS OF NODE LINES, 


EDGE CONDITIONS, AND MATERIAL CONSTANTS 
(a) Invar panel 


Location of node lines 

Measured natural frequencies, Hz 

Theoretical natural frequencies for 
edge restraints, q^ = 14, q^ = 110.0, Hz 

Deviation of measured frequencies from 
theoretical values, percent 

Theoretical natural frequencies for 
clamped edges, q^ = q^ = °°, Hz 


Deviation of measured frequencies from 
theoretical values, percent 


18.1 

Panel characteristics 


38.0 

38.9 

-2.3 

46.4 






1 i 1 

1 i i 1 


liiil 

52.0 

79.0 

106.0 

116.0 

53.1 

77.9 

102.9 

112.9 

-2.1 

1.4 

3.0 

2.8 

60.1 

84.4 

120.9 

119.6 

13.5 

6.4 

12.3 

3.0 


Handbook 


Panel dimensions: 9.00x18.00x0.0192 in. 

Edge conditions: q x = 14, q y = 110.0 

Invar, 64 percent Fe, 36 percent Ni 
21.0xl0 6 lb/in. 2 
0.29 

507.9 lb/ft 3 
0.7xl0' 6 in./in./°F 

(b) Magnesium panel 
Location of node lines 


Material 
E: 
v : 

Ps : 

a : 


Measured natural frequencies, Hz 

Theoretical natural frequencies for 
edge restraints, q^ = 62, q^ = °°, Hz 

Deviation of measured frequencies from 
theoretical values, percent 

Theoretical natural frequencies for 
clamped edges, q^ = q^ = Hz 

Deviation of measured frequencies from 
theoretical values, percent 

Panel characteristics 



1 . 1 






1 

1 

LLU 


1 i i ! 1 

110.0 

143.0 

212.0 

287.0 

298.0 

110.1 

144.9 

207.0 

287.7 

295.5 

-.1 

-1.3 

2.4 

-.2 

.9 

115.7 

149.8 

210.3 

301.2 

297.9 

-4.9 

-4.5 

.8 

-4.7 

.03 


Panel dimensions: 
Edge conditions: 
Material : 

E: 

v: 

p s : 

a : 


9.00x18.00x0. 0401 in 
q x = 62, q = « 

Magnesium alloy AZ 31B-H24 
6.2xl0 6 lb/in. 2 ' 

0.35 

110.0 lb/ft 3 
14.2X10 -6 in./in./°F 


Measured, see p. 3 
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TABLE 2.- TABULATED PANEL FLUTTER DATA INVAR PANEL 


Mach 

number 

Wall 

displacement 

Free-stream 

temperature 

Dynamic pressure for 
multiples of 
subflutter response, 
psf 

Boundary- layer 
thickness, in. 

Flutter 
frequency , 
Hz 

Reduced 

flutter 

frequency 

Supersonic 

reduced 

flutter 

frequency 

Mass 

ratio 

parameter 

Cavity 

density 

ratio 

Modified 

stiffness 

parameter 

Stiffness 

parameter 


H, in. 

T ro , °R 

^2X 

q 

n 3X 



q 6X =q F 

6 

6 m=i.o 

6* 

0 

*F 

k 

K 

y 


k* (1) 

k ( 2 ) 

1.20 

I 

0 

— 

416.4 

370 

382 

392 

400 

403 

1.14 

0.52 

0.228 

0.110 

37 

0.145 

0.396 

0.0166 

0.0133 

0.0309 

0.0252 

1.20 

.30 

412.5 

273 

304 

315 

321 

325 

.82 

.32 

.162 

.076 

42 

.166 

.452 

.0135 

.0107 

.0310 

I .0253 

; 1.20 

i 

.50 

415.2 

243 

(251 254 

255 

256 

.66 

.26 

.134 

.058 

45 

.177 

.483 

.0106 

.0084 

1 .0309 

! .0252 

! 

1.20 

| .80 

411.4 

208 

215 

217 

219 

221 

.46 

.25 

.108 

.045 

46 

.182 

.496 

.0092 

.0073 

t 

1 .0311 

1 .0254 

, 1.30 

! o 

] 396.8 

423 

•443 

451 

458 i 464 

1.15 

.54 t .262 

.112 

47 

.175 

.329 

,0171 

j .0131 

1 .0316 

.0258 

, 1.30 

.80 

397.0 

267 

,275 

277 

279 f 280 

.38 

.16 

.102 

( .038 

52 

' .193 

,364 

.0103 

j .0078 

.0316 

.0258 

! 1.30 

1.00 

396.6 

1 

261 1 263 

264 

265 

265 

.32 

,14 

.088 

.0381 53 

! 

1 .197 

1 

.371 

.0098 

j .0074 

1 .0316 

.0258 

1.40 

I .60 

i 381.3 

1373 

'388 

390 

t 

'390 

390 

* .47 

.12 

.105 

.042 

48 

| .169 

.246 

j .0129 

1 

.0095 

| 

.0323 

.0263 

j 

1.40 

.80 

381.0 

349 

'366 

I370 



|370 

370 

! .34 
.1 

.12 

.094 

) 

.033 

i 

52 

.183 

.267 

1 

.0123 

I 

1 

.0090 

1 . _ 

.0323 

J 

.0263 


: E = 21.0xl0 6 lb/ in. 2 (Handbook) 
2 E = 14.1xl0 6 lb/in. 2 (See p. 8) 


Hr- 0.120 


0.020 


Note: All dimensions 



Figure 1.- Invar panel 


H - -^1 



and construction details 








(a) Three -fourths front view of airstream side 
Figure 3«- Variable boundary-layer test fixture. 


A-39842 


NJ 

4 ^ 


(b) Three -fourths rear view of plenum chamber side 


A-39841 



Figure 3*- Concluded. 



to 

in 


Figure 4 .- Panel with variable houndary-layer fixture installed in the 2-hy 2-foot 

transonic wind tunnel. 


A-35851. 1 



Figure 5*“ Typical variation of maximum panel response with dynamic pressure 





Figure 7.- Dynamic pressure for multiples of subflutter response as a 
function of boundary- layer thickness. 



\0\tQ 



Figure 8.- Normalized flutter dynamic pressure as a function of normalized 

boundary- layer thickness. 





5=0 


Extrapolation of linear 
approximation of flutter 
frequency ratio 


• 2h 


0 *■ 

0 


.02 .Oh .06 


a 


Figure 10.- Normalized flutter frequency as a function 



the normalized boundary- layer thickness 







0 .02 .04 .06 .08 .10 .12 
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(a) Comparison at = 1.30 

Figure l4.- Comparison of experimental flutter points with theoretical 
stability boundaries in the stiffness parameter -mass ratio (k, \i) 
plai 





Theoretical flutter 
boundaries 
g a/b 

0I05 o '.46 I Ref ‘ 2 ^ mode anal y sis ) 



(c) Comparison at = 1.20 


Figure 14 . - Continued . 
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